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We study the spin-lattice relaxation rate of nuclear magnetic resonance in a two-band supercon- 
ductor. Both conventional and unconventional pairing symmetries for an arbitrary band structure in 
the clean limit are considered. The importance of the inter-band interference effects is emphasized. 
The calculations in the conventional case with two isotropic gaps are performed using a two-band 
generalization of Eliashberg theory. 
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I. INTRODUCTION 



Although the Fermi surface in most superconductors 
consists of more than one sheet, this does not neces- 
sarily mean that all those materials are multi-band su- 
perconductors. The true multi-band (in particular two- 
band) superconductivity is in fact a rather uncommon 
phenomenon characterized by a significant difference in 
the order parameter magnitudes in different bands. For 
this to be the case, the system has to satisfy some quite 
stringent requirements, namely the pairing interactions 
and/or the densities of states should vary considerably 
between the bands and the inter-band processes, e.g. due 
to impurity scattering, should be weak. Although some 
examples have been known since early 1980s^ the recent 
swell of interest in this subject has been largely stimu- 
lated by the discovery of two-band superconductivity in 
MgB2. 2 Most of the experimental evidence, see Ref. 3 
and the references therein, support the conclusion that 
there are two distinct superconducting gaps A n and A^ 
in this material, with A^/A^ ~ 2. 63. 4 (There are actu- 
ally four bands crossing the Fermi level in MgB2, which 
can be grouped into 2 quasi-two-dimensional cr-bands and 
2 three-dimensional 7r-bands and described by an effec- 
tive two-band model.) Other candidates for multi-band 
superconductivity which have emerged recently include 
nickel borocarbides^ NbSe2i^ and also the heavy- fermion 
compounds CeCoIn 5 (Ref. 0) and CePt 3 Si (Ref. i). It 
seems more likely to find a two-band superconductivity in 
unconventional materials, since they are intrinsically in a 
clean limit, so at least the gap averaging due to impurity 
scattering is not effective. 

Theoretically, a two-band generalization of the 
Bardeed-Cooper-Schrieffcr (BCS) model was intro- 
duced independenly by Suhl, Matthias, Walker^ and 
Moskalenko. 10 In subsequent developments, many as- 
pects of the multi-band model have been studied, in- 
cluding the thermodynamic and transport properties, the 
effects of impurities and strong coupling, efc>ii*i2*i2ii4 
Surprisingly, little attention has focussed on such an im- 
portant characteristic as the spin-lattice relaxation rate 
of nuclear magnetic resonance (NMR). The mea- 
surements of T^ 1 probe the properties of the electron 
subsystem which are local in real space and therefore ex- 



tremely non-local in the momentum spaced In the pres- 
ence of multiple Fermi-surface sheets this would give rise 
to inter-band interference terms in Tj~ , even without 
any inter-band scattering due to interactions or impuri- 
ties. The inter-band terms in T^ 1 are not negligible and 
can be expected to strongly affect the temperature depen- 
dence of the relaxation rate compared to the single-band 
casei 1 ^ 

The purpose of this article is to calculate the nuclear 
spin relaxation rate in a two-band superconductor, for 
both conventional and unconventional types of pairing. 
We focus on singlet pairing in the absence of impuri- 
ties, assuming that the relaxation is dominated by the 
Fermi contact interaction between the nucleus and the 
conduction electrons. The article is organized as follows. 
In Sec. [H] we develop a general formalism based on an 
anisotropic two-band BCS model and show that, while 
the resulting expressions in the unconventional case are 
well-defined and can be calculated without any additional 
complications, in the conventional isotropic case one en- 
counters divergent integrals. In Sec. 11111 we single out the 
isotropic case for a strong-coupling theory treatment, in 
which the divergences are smeared out due to the quasi- 
particlc lifetime effects. In Sec. IIVI we apply the general 
theory to the relaxation rate on the 25 Mg site in MgB2 
using the realistic strong-coupling parameters. 



II. WEAK COUPLING THEORY 

Without the loss of generality, we consider the case of 
a nuclear spin 1 — 1/2 located at the origin of the crystal 
lattice. Higher values of / change only the overall pre- 
factor in the expression for the relaxation rate^ which 
drops out of the ratio of the relaxation rates in the su- 
perconducting and the normal states. The spin-lattice 
relaxation rate due to the hyperfine contact interaction 
of the nucleus with the band electrons is given by 
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R=^==- J - lim *»*f-M, 
T\T 2ir wo— >o ujq 



(1) 



where J is the hyperfine coupling constant, ujq is the 
NMR frequency, and K^_(u>o) is the Fourier transform 
of the retarded correlator of the electron spin densities 
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at the nuclear site: 



K*{t) = -i([S + (0,t),S-(tO,Q)])W)- 



(2) 



Here S±(r,t) = e 
Hamiltonian, and 



iH.t i 



S±(r)e 1 e , H e is the electron 



5 + (r) = ^|(r)Vi(r), S-(r) = i>\(r)^(r) (3) 

(h = kg = 1 in our units, and the spin quantization 
axis is along the external magnetic field H) . The deriva- 
tion of Eq. (JTJ is outlined in Appendix The re- 
tarded correlator is obtained by analytical continuation 
of the Matsubara time-ordered correlator: K+_(ui) = 
K(v m )\ ll/m _ UJo+i0 +, with v m = lirmT. 

We assume that there are two spin-degenerate elec- 
tron bands in the crystal (the generalization to an arbi- 
trary number of bands is straightforward), and neglect 
the spin-orbit coupling. The two-band generalization of 
the BCS Hamiltonian reads H e 



H + H int jL where 



Hn = 



E 

i.ka 



(4) 



breaking (p a 's can be chosen real. The basis functions do 
not have to be the same in both bands, but we neglect 
this complication here. 

The properties of our superconductor can be described 
using a standard field-theoretical formalism in terms of 
the normal and anomalous Gor'kov functions^ 

G i>a p{k, t) = S a pGi(k, r), 

Fl ap (k,T) = {-ia 2 ) aP F}{k,T), 

which can be combined into a 2 x 2 matrix Green's func- 
tion 

Gi(k ' T) -[-FHk,T) - Gi (-k,-T))- (?) 

In the mean-field approximation, the interaction Hamil- 
tonian is reduced to the form 



Hint = \ E A ^4 M c[_ kl + H.c. 



(8) 



i,k 



is the non- interacting part [i = 1,2 is the band index, 
a =|, { is the spin projection, and the chemical potential 
\x is included in the band dispersion) , and Flint = H. 



(i) 



H {2 > + H 



(12) 



is the pairing interaction. For anisotropic 



singlet pairing, we have 

k,k' 

H int =\^V 12 {k,k l )c\ M c\_ ki c 2 ,-k>iC2,k'] 

+H.c. (5) 



fe.fe' 



The Hamiltonians and Fc^l describe the intra-band 

(12) 

pairing of electrons, while Fc\ nt describes the pair scat- 
tering between the bands. The inter-band interactions 
of the form c\ k ^c\ _p l ^C2,-k'lCi,k'l are suppressed if the 
band splitting is large compared to all energy scales rel- 
evant to superconductivity. We assume, following Hebel 
and Slichter— that, while the resonance is observed in 
a strong field in the normal state, the relaxation takes 
place in a uniform superconducting state after switching 
off the field. 

The pairing symmetry is the same in both bands and is 
determined by one of the irreducible representations, T, 
of the point group of the crystal. The functions Vy (k,k') 
are non-zero only in a thin energy shell near the Fermi 
surfaces and can be represented in a factorized form: 



dr 



V^ti) = V ij Y, t Pa{k)y a (k'), 



(6) 



a=l 



where tp a (k) are the basis functions, and dp is the dimen- 
sionality of r. In the absence of time-reversal symmetry 



where A^^ is the superconducting order parameter in the 
ith band, which can be written as 



Ai jfe = y^7?i,ay a (fc), 



(9) 



with rji ta being the order parameter components. Both 
order parameters appear at the same critical temperature 
T c , but have different temperature dependences, which 
can be found by solving a system of 2dr self-consistency 
equations for the functions rji ia (T). In the frequency rep- 
resentation, the Green's functions (0 become 

Gi ^ = - .g+eL + lA^l 2 ' (10) 

where r, are Pauli matrices, w„ = (2n + 1)tvT, and 



X. k = [ £ k *J fc ]■ ill) 



Now we return to the calculation of the relaxation rate 
(|T|). For zero spin-orbit coupling, the spin operators © 
can be written in the band representation, using 



v i.k 



(12) 



where Ui.k{r) are the Bloch functions, which are pe- 
riodic in the unit cell, and V is the system volume. 
Inserting these into Eqs. J3J), one obtains the Mat- 
subara spin correlator K(t) = —(F t S + (O,t)S-(O,0)), 
which can be decoupled in the mean-field approxima- 
tion, using the Green's functions Q. In the absence 
of time-reversal symmetry breaking, one can show that 
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Ui t —k(0) = u* k {0). Then, taking the thermodynamic 
limit, we have 



Un + V m )G{k 2 ,UJn)}, (13) 



where 



and 



d D k 



/(-) -Ji 7 ?<->-/£ 



g{k,u n ) = ^K fe (o)| 2 G i (fe,( 



(•••), 



(14) 



with Gi(k,u> n ) given by Eq. (|10fl . 

Calculating the matrix traces and the Matsubara sums 
in Eq. lfH3|) followed by the analytical continuation to real 
frequencies, we find that the imaginary part of K+_ (u) 
is proportional to to at u> — > 0. The momentum integrals 
are calculated making the usual assumption that u^fc(O) 
and A^fe weakly depend on in the vicinity of the 
Fermi surface (i.e. within the energy range of the order 
of T). We introduce the local density of quasiparticle 
states at r = 0: N(u) = Ni(w) + N 2 {uj) (u > 0), where 



Ni{u) = - I \u iik {0)\ z 8(uj - E iik ) 



K,fe(0) 



|Ai, fe | 



> (15) 



where Ei >k = y £ 2 fc + \&iM 2 ^ s tne Bogoliubov excita- 
tion energy in the ith band, the angular brackets stand 
for the average over the Fermi surface, and Np,i — 
(1/87T 3 ) J dSp/\vF,i\ is the density of states at the Fermi 
level in the ith band. The angular integration in Eq. 
((TKjl is restricted by the condition |A^fc| < uj. We also 
introduce the function M(uj) = Mx{u) + M^ui), where 



i r Ai 



2 J Ei.k 

k 



\u hk (0)\ 2 6(cu-E hk ) 



N 



FA 



K,fe(0) 



H.k 



V^ 2 " lAi.fcl 1 



(16) 



Then, 



Of 



R = J 2 Jdoj \-ZL) [N 2 (oj) + \M(cu)\ 2 ] , (17) 
o 

where f(u>) — (e"/ T + is the Fermi function. 

For Aj(fe) = 0, we have M(u>) — 0, and the normal- 
state relaxation rate is given by R n — J 2 N 2 /2, where 



N n>i =N Fii (\u i>k (0)\ 2 ) 



(18) 



Finally, we obtain for the ratio of the NMR relaxation 
rates in the superconducting and the normal states 



R a n 
flu 



oui 



N 2 (co) + \M(cj)1 



(19) 



As we pointed out at the beginning of this Section, our 
result does not depend on the nuclear spin /. The expres- 
sion Ijl9|l has two notable properties. First, the relaxation 
rate is controlled by the local densities of quasiparticle 
states. Only in the limit of a single-band isotropic pair- 
ing can one express R in terms of the total density of 
states and recover the Hebel-Slichter formula^ see Sec. 
Ill AI below. Second, the contributions to the spin-lattice 
relaxation rate from different bands are not simply addi- 
tive, since there are inter-band interference terms in Eq. 
()19J) . These terms are present even in the absence of any 
inter-band interactions or impurity scattering and can be 
traced back to the local character of the hyperfine cou- 
pling IS, which mixes together the electron states near 
the Fermi surface from different bands. 



A. Conventional pairing 

The order parameter is "conventional" if it trans- 
forms according to the unity representation of the point 
group The gap functions A^*. can be isotropic or 
anisotropic, with M s (oj) ^ in both cases. 

Assuming the isotropic pairing with a uniform order 
parameter, we have A^j, = Aj, where both gap functions 
can be chosen real without loss of generality. One can 
view this as an extreme case of anisotropic superconduc- 
tivity on an extended single sheet of the Fermi surface, in 
which the gap function is allowed to take only two values, 
Ai and A2. The densities of states become 



N(u) 



N„ 



M(w) = J2N n 



(20) 



(21) 



Substituting these expressions in Eq. (|19fl we arrive at 
a logarithmically divergent integral. The origin of this 
divergence is the same as in the Hebel-Slichter formula 
in the single-band case^ one has to square the BCS- 
like density of quasiparticle states, which is singular at 
E = Ai, A2. Allowing for a non-zero NMR frequency Uq 
yields the relaxation rate which is still much higher than 
that observed in experiment^ 

One can smear out the singularity and cut off the di- 
vergence either by introducing some gap anisotropy^ or 
by taking into account the strong-coupling effects, which 
lead to a finite lifetime of quasiparticles and therefore to 
energy-dependent complex gap functions. 21 Which mech- 
anism is more important depends on the material. In Sec. 
nil below, we adopt the latter point of view and derive 
the strong-coupling expression for the relaxation rate for 
an isotropic gap. 
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B. Unconventional pairing 

If the order parameter transforms according to a non- 
unity representation of the point group, then it fol- 
lows from the obvious property of the Bloch functions 
\ui, g k(0)\ 2 — |u ii fc(0)| 2 (g is an arbitrary element of the 
point group) that M(u>) = 0. Therefore, 



|i = 2 / duj 

Rn J 



dui 



Nijuf) + N 2 (uj) 
N n<1 + N n . 2 



(22) 



where Ni(w) and N n< i are defined by Eqs. (|15|l and l|18|) 
respectively. In most cases the integral converges, be- 
cause the square-root singularity in the density of states 
is smeared out by the intrinsic anisotropy of the gap. 
The only exception is an unconventional order parame- 
ter with isotropic gap (e.g. an analog of the B-phase of 
3 He in a charged isotropic superfluid), in which case the 
integral is again logarithmically divergent. 

Since the inter-band pair scattering in Eq. (JSJ 

induces the order parameters of the same symmetry in 
both bands, the low-energy behavior of Aq(w) and N 2 (uj) 
is characterized by the same power law. If there are line 
(point) nodes in the gap, then Ni >2 (u>) oc w (u> 2 ) at uj — > 
0^ and R oc T 2 (T 4 ) as T — > oM^. This behavior has 
indeed been observed in most heavy- fermion compounds, 
for a recent review see Ref. I24I 

This picture will change if the gap magnitudes in the 
bands are considerably different (as mentioned in the In- 
troduction, there are indications that this might be the 
case in such materials as CeCoIns and CePtaSi). For ex- 
ample, if the gap in one band is much smaller than in 
the other, then, taking the limit A 2t k — > 0, one obtains 
instead of Eq. J22J 



Rn 



N 2 (cu) + 2N 1 {u;)N n , 2 + N 2 
(JV n ,i + N n , 2 ) 2 



(23) 



While the last term in the integral contributes to the 
residual relaxation rate at T = 0, it is the second term 
that controls the power-law behavior at low T: we now 
have R = const + aT for line nodes, and R = const + aT 2 
for point nodes. 

As an illustration of the above results, let us consider 
a simple example of a quasi-two-dimensional two-band 
superconductor with circular Fermi surfaces and a d-wave 
gap Ai,fc = Aocos26', which has vertical lines of nodes. 
The fraction of the density of states from the electrons in 
the unpaired band is r — N Ui2 / (N nt i+N n ^ 2 ). The Fermi- 
surface average in Eq. 115|) can be done analytically: 



xK{x 2 ) 
K(x- 2 ) 



if x < 1, 
if x > 1, 



(24) 



where x — w/Aq, and K(x) is the complete elliptic inte- 
gral of the first kind^ 




FIG. 1: The NMR relaxation rate in a two-band supercon- 
ductor with lines of nodes in one band and a negligible gap 
in the other, for different values of r = N n , 2 /(N n ^ + AT„ i2 ). 



In Fig. ^ we show the results of the numerical cal- 
culation of the temperature dependence of the relax- 
ation rate (|23|l for different values of r. Instead of de- 
termining the exact temperature dependence of Ao at 
all T, which would involve a full numerical solution of 
the self-consistency gap equation, we use the a pproxi- 
mate expression A o (T)/A o (0) = y/l - {T/T c ) 3 , where 
Ao(0)//cbT c = 1.30 (this number is obtained from the 
solution of the gap equation at T = 0). For r = 0, one 
recovers the limit of a single-band d-wave superconductor 
with R oc T 2 at low T and a small Hebel-Slichter peak 
immediately below T c . As r grows, so do both the devi- 
ation from the T 2 behavior and the residual relaxation 
rate at T — 0. One interesting observation is that even 
if the density of states is dominated by the contribution 
from the unpaired sheet of the Fermi surface, one still 
can see an appreciable suppression of the relaxation rate 
at low temperatures. 



III. STRONG COUPLING THEORY 

In this section we generalize the results of the weak 
coupling theory, Sec. II, to the case of an electron- 
phonon multi-band superconductor which could be de- 
scribed by Eliashberg-type equations* 4 ^ To include the 
self-energy effects associated with both electron-phonon 
and screened Coulomb interaction one replaces Eq. I|13(l 
with 

K(v m ) = ±T y £ J Eki, fel (0)| 2 |% )fe2 (0)| 2 

x Tr [d (fci , w n )Tij (fci , fe 2 ; 0J n , v m ) 

xGj(k 2 ,Un + Vm)], (25) 
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where Gi(k,u) n ) are given by 

Gi(fc ' w " j - <#?>») + tf lfc + ^ fc K) ' 

(26) 

instead of Eq. (|10|l . Here Zi ik (u> n ) and (j>i,k(un) are 
the renormalization function and the pairing self-energy, 
respectively, for the «th band. 

The vertex functions Tij(ki,k2',oJ n ,^m) = 
Tij{k,uj n ;q 1 v m ) need to be calculated in the conserving 
approximation consistent with the approximations 
used to calculate the electron self-energiesi 26 i 27 i 28 Since 
after analytic continuation iv m — > luq + i0 + one is 
interested in the low- frequency limit, see Eq. l|T]l. 
and the Migdal's theoremSLSi guarantees that the 
electron-phonon contribution to the vertex functions 

satisfies lim^^o f \j ^ {k, u n \ <Z, v m ) ~ tq for any finite 
q, the electron-phonon interaction can be suppressed in 
evaluating the vertex parts. The Coulomb interaction, 
on the other hand, leads to Stoner-type enhancement^ 
which is unaffected by the superconducting transition 
(assuming the usual electron-phonon pairing mecha- 
nism) and thus should cancel out from the ratio R s / R n . 



Hence, we replace f y in Eq. I|25[l with the unit matrix 
To in computing the ratio of the spin-lattice relaxation 
rates in the superconducting and normal states. We 
note, however, that the single particle energies & & are 
assumed to be renormalized by the Coulomb interaction 
and that the electron-phonon vertices entering various 
self-energy parts in Gi(k,Lu n ) are Coulomb vertex 
corrected and Coulomb renormalized as discussed in 
Ref. |H 

Next, one introduces the spectral representation for 
Gi(k,uj n ) 

+00 

G i (fe, w „)= / <ko$^-, (27) 

J lLU n -UJ 
— 00 

with 

Ai(k,uj) = --lmG l (k,uj + iO + ), (28) 

7T 

which allows one to calculate the Matsubara sums in Eq. 
(|25|l . followed by the analytical continuation iv m — > luq + 
i0 + . In the limit luq —> we obtain 



1 ImKfujQ - 
hm 

Wo— >Q 7T UJQ 



i0 + ) _ 1 



Im 



Im 



7T 

fci.a — 00 

ujZ jM (u) 



Im 



^(-£)Ei^(°)i 2 i u W°) 

i,3 



L0Zi. kl {ui) 

A,fei(w) Dj. k2 {uj) Di. kl {uj) D jtk2 {u) D iM {uj) D jtk2 (u!) 



Im ■ 



Im 



(29) 



where 



Di. k {u) = [uZ ith (u)Y 



(30) 



and Z iik (u) = Zi^ k {uj + i0 + ), <f>i,h{w) = &,fe(w + i0+). 

Next, we assume that Zi }k (uj) and 0i,fc(w) are isotropic, 
which seems to be a reasonable assumption for MgB2, 4 
and use a weak dependence of these functions on ^ ^ 
which is one of the consequences of the Migdal's theorem. 
Hence, the fe-dependence of Zi and & can be suppressed, 
and after defining the local densities of states (|15|l . (|16|l . 
and l|18(l , the momentum integrations in Eq. 129(1 can be 
easily performed. The final result has the form 



+00 

^ = 2 fduf-lf) ^)+f 2 n 

Rn J v Nl 





where 



JV(w) = ^A^Re 

i 

M(lo) = J2N„,iRe 



Ai(w) 



\A> 2 -A 2 H 



(31) 

(32) 
(33) 



and Ai(w) = <pi(u) / Zi{uj) is the gap function in band 
i. In arriving at l|31|) we have used A^(— ui + i0 + ) = 
A*(cj + z0 + ) which follows directly from the spectral rep- 
resentation (|27|) . It is easy to sec that our result l|31|l . 
(|32|l . i|33|) reduces to the one given by FibichSi in the 
case of a single isotropic band, and to Eqs. I|19fl . 120fl . 
(|21|l in the weak coupling limit, when the gap function 
does not depend on lj. Similar to the single-band case, 
the presence of non-zero imaginary parts in Ai(w) leads 
to the smearing out of the BCS square-root singularities 
in JV(w) and M(u>). 



IV. APPLICATION TO MgB 2 

For a quantitative application of the results of the pre- 
vious section to a particular compound, one needs to 
know both the band-structure characteristics and the in- 
teraction parameters of the Eliashberg theory. The only 
two-band superconductor for which these are presently 
available is MgB2. 

Different contributions to the hyperfine interaction in 
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Mgl$2 were calculated using the local-density approxi- 
mation in Refs. I33ll34i It was found that, while the re- 
laxation at the 25 Mg nucleus is dominated by the Fermi 
contact interaction, for the 11 B nucleus it is the inter- 
action with the orbital part of the hypcrfinc field that 
makes the biggest contribution. These predictions were 
subsequently found to be in excellent agreement with ex- 
periments in the normal state. 35 - 36 - 37 To the best of our 
knowledge, the experimental results on temperature de- 
pendence of Tj -1 in the superconducting state of MgB 2 
are available only for the 11 B nucleus ^^JSdL There- 
fore our theory, which should be applicable only to the 
relaxation rate for the 25 Mg nucleus in a clean sample, 



where 



With a set of four electron-phonon coupling functions 
a 2 Fij(Q), i,j — a, n, calculated in Ref. 0, and with a set 
of the Coulomb repulsion parameters fj,*j(u> c ) determined 
in Ref. 42 to fit the experimental critical temperature T c , 
Eqs. 1341351) were solved for the complex gap functions 
A cr (w) and A„(uj) at a series of temperatures below T c . 
A representative solution near T c is shown in Fig. [2] (T = 
0.968T C ). 

The band structure calculations 43 indicate that the 
contribution to the local density of states at the Mg site 
from the a band is much smaller than that from the 7r 



cannot be directly verified by comparison with the exist- 
ing experimental data. The lack of data on Tj -1 for the 
25 Mg nucleus is presumably related to the small magnetic 
moment and a low natural abundance of this nucleus as 
discussed in Re f. l35i Nevertheless, the experiments per- 
formed in Refs. ISFillSfil indicate that it is possible in princi- 
ple to measure 25 R below the superconducting transition 
temperature. 

To calculate R s /R n in the superconducting state of 
MgB2 we have solved the coupled Eliashberg equations 
with the realistic interaction parameters for the isotropic 
two-band model^ on the real frequency axis and at finite 
temperature: 



(34) 



(35) 



(36) 



(37) 

I 

band. Therefore we can set N a = in the expressions for 
T^ 1 on the 25 Mg nucleus. In Fig.[3]we show the temper- 
ature dependence of R s / R n obtained from the numeri- 
cal solutions of the strong-coupling gap equations, using 
Eqs. (|31l I32I33|) . At the lowest temperatures, the relax- 
ation rate is exponentially small, while at T — ► T c — 0, 
R s /R n - 1 is proportional to (1 - T/T c )°- 5 . The most 
prominent qualitative feature is a shift of the Hebel- 
Slicher peak away from T c to a lower temperature, at 
which the coherence factor from the lower gap in the tt- 
band makes the maximum contribution. The significant 
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A,- (a/) 



a/ 2 -A 2 (a/) 



f{-ur>K+ M (u,tJ) - f(u>')K +tij (u,-u/) 



(cj c ) tanh — + Kl% (u, u/) - K$% (w, -u/) 



+ 00 



Zi(w) = l-iV / dJ Re [f(-uj')K_ tij (uj, u/) ~ f(u')K- Aj (u, -</) 
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FIG. 2: The solutions for the real and imaginary parts of 
Aa-(w) and Att(uj) in the entire phonon energy range for 
MgB 2 , at T = 0.968T C . The inset shows the solutions in 
the low energy range where the real parts of the gaps are 
quadratic functions of uj and the imaginary parts of the gaps 
are linear functions of uj at low enough energy for T >0. 




T/T 



FIG. 3: The ratio R„/R„ as a function of the reduced tem- 
perature T/T c in the case when the relaxation is dominated 
by the lower-gap band. 



nuclear spins is dominated by the Fermi contact inter- 
action with the band electrons. Our main result is that 
there are important inter-band contributions not related 
to any scattering processes, which change the tempera- 
ture dependence of the relaxation rate. In particular, if 
there are unpaired sheets of the Fermi surface in a su- 
perconductor with gap nodes, then in addition to the 
residual relaxation rate at T = 0, one should see unusual 
exponents in the power-law behavior at low T. The ob- 
servation of those exponents could be a strong argument 
in favor of multi-band superconductivity. 

To illustrate the general theory, we calculated the re- 
laxation rates in the clean limit for (i) a two-dimensional 
d-wave superconductor, using the BCS theory, and (ii) 
an isotropic s-wave superconductor, for which a strong- 
coupling treatment is required. In the latter case, we ap- 
plied our model to the 25 Mg nucleus in MgB2, for which 
the relaxation is due to the Fermi contact interaction 
and the parameters of the Eliashberg theory are known. 
The predicted temperature dependence of the relaxation 
rate is quite unusual and should be easily detectable in 
experiments. 

In order to expand the applicability of our theory, one 
should include disorder, especially the interband scat- 
tering, which is a pair-breaker in the multi-band super- 
conductors. Although the unconventional candidates for 
multi-band superconductivity, such as CeCoIns, are in 
the clean limit, in general the impurity effects might be 
significant. Also, our basic assumption that the relax- 
ation is controlled by the local fluctuations of the Fermi- 
contact hyperfine field, can be violated in some cases, e.g. 
for the 11 B nucleus in MgB2. Another possible general- 
ization would include the effects of the gap anisotropy 
within the separate bands4& It is well knowni^ that the 
spread in gaps within a single band leads to the suppres- 
sion of the coherence peak in R s /R n below T c . Finally, if 
the NMR measurements are done at a non-zero magnetic 
field in the presence of vortices, then the inhomogencity 
in the order parameter in the mixed state strongly af- 
fects the density of quasiparticle states and therefore the 
relaxation rate^I 
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V. CONCLUSIONS 

We calculated the NMR relaxation rate T^ 1 in a sin- 
glet two-band superconductor without spin-orbit cou- 
pling and impurities, assuming that the relaxation of the 



APPENDIX A: DERIVATION OF EQ. © 

We assume that the dominant mechanism of the spin- 
lattice relaxation is the interaction between the nuclear 
spin magnetic moment fvy n I (7„ is the nuclear gyromag- 
netic ratio) and the hyperfine field created at the nucleus 



8 



by the conduction electrons. The system Hamiltonian is 
H = H e + H n + H int , where H e describes the electron 
subsystem, H n — —fvy n IH is the Zeeman coupling of the 
nuclear spin with the external field H , and 



H:, 



-hj n Ih 



(Al) 



is the hyperfine interaction. For / = 1/2, we have two 
nuclear spin states I z = ±1/2 with the energies Ej, = 
—fkuoI Zl where ujq = j n H is the NMR frequency and the 
spin quantization axis is chosen along H . The hyperfine 
field h can be represented as a sum of the Fermi contact, 
the orbital, and the spin-dipolar contributions 15 . Their 
relative importance depends on the electronic structure 
and therefore varies for different systems. For example, 
if the Fermi contact interaction is dominant, then h = 
— (8n/3)h-f e S(0), where 7 e is the electron gyromagnetic 
ratio and S(r) = (l/2)cr a pTpl i (r)Tpp(r) is the electron 
spin density at r = 0. The derivation below does not 
rely on any particular expression for the hyperfine field. 

According to Ref. fl5L the relaxation rate for a spin- 1/2 
nucleus is given by 



1 



= W+- + W— 



(A2) 



where W-\ and W |_ are the transition probabilities 

per unit time, from I z = +1/2 to I z = —1/2 and from 
I z = —1/2 to I z = +1/2, respectively. The hyperfine in- 
teraction is usually small, which makes it possible to use 

the lowest-order perturbation theory to calculate W-\ 

and W-+. The states of the whole system at zero hyper- 
fine coupling can be represented as |J) = \i,I z ), where 
i labels the exact (in general, many-particle) eigenstates 
of H e , with energies E^. When J ^ 0, then the transi- 
tion probability per unit time from an initial state \I) of 
energy £j to a final state \F) of energy £p can be found 
using the Golden Rule: 

w\i)^\F) = Y|(/|/f int |F)| 2 5(^ - £ F ). (A3) 

The transition rates for the nuclear spin are calculated 
in the usual fashion by averaging over the initial and 
summing over the final electron states. 

For VF + _, we have \I) = \i,+l/2), £/ = E, - tkJ /2 
and \F) = |/, -1/2), £ F = E f + fiw /2. Then 

W+- = ^ Pe,i U>|<,+l/2>-|/,-l/2> > ( A4 ) 
/ 

where p e — e~P H ' /Tr e~~P H " is the density matrix of the 
electron subsystem. Inserting here the expressions (IA3|) 
and (|A1() and representing Ih = I z h z + (I + h_+I_h + )/2, 



where I± = I x ±il y and h± — h x ±ih y , we find that only 
the I-i-h_ term makes a non-zero contribution. Therefore, 



W+- = 



YPe,i\{i\h-\f)\ 2 S{Ei-E f -fuj 



This expression can be simplified by using the identity 



f°° dt 

5(E i -E f -n, o) = ]_ oa —< 



i(Ei-Ef-hui )t/h 



and the fact that h_ — h + , which allow us to write 

\(i\h_\f)\ 2 e^ Ei - E fW n 

= (i\e iE ^ h h^e- iE f^ h \f)(f\h + \i) 

= {i\h-(t)\f){f\h + (0)\i), 

where h±(t) = e lH ^ t / h h ± e~ lH " t ^ h . Now the sum over the 
final states can be calculated, and we finally have 



W +- = f 



dt e 



-iujQt 



(M*)Mo)>. 



(A5) 



The angular brackets here stand for the thermal averag- 
ing with respect to the electron density matrix p e . Simi- 
larly, we obtain 



W_+ = ^/ dte lwot (h+{t)h^(0)). (A6) 



Combining Eqs. (|A5|) and JAB}, we have 

k = 1 f /^* e " ot({Mt) ' M0)}) - (A7) 

The integral on the right-hand side here can be ex- 
pressed in terms of the Fourier transform of the re- 
tarded correlator of the hyperfine fields K^ h (t) — 
-i{[h+(t),h-(Q)])0(t), giving 



ll k B T 
2ir hu>o 



(A8) 



Here we used the fact that in a typical experiment the 
condition Huq <C feT is always satisfied [we also note 
that since W+-/W-+ = e - 0huJO ~ 1 due to the detailed 
balance in the thermal equilibrium, one could use T^ 1 ~ 

2W1| instead of l|A2|l ]. Keeping only the Fermi contact 

term in the hyperfine interaction (|A1|) . we finally arrive 
at Eq. P with J = (8n/3)j n j e . 
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